Relating monomer to center-of-mass distribution functions in macromolecular fluids: 

extension to arbitrary systems 
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In the framework of the polymer reference interaction site model (PRISM) extended to incorporate 
the centers of mass (CM) of the groups of monomers constitutive of the macromolecules as auxiliary 
sites, a general relationship between the monomer-monomer and CM-CM pair correlation functions 
is derived for arbitrary mixtures of macromolecules of arbitrary structure. It extends previous results 
for homopolymer fluids and should be useful to develop semi-analytic coarse-grained descriptions of 
complex fluids. 



Following the recent successes of a mesoscale coarse- 
graining strategy for polymeric systems, in which the 
microscopic chain configurations are mapped onto those 
of a fluid of mesoscopic particles by the introduction of 
effective state-dependent pair interactions between the 
centers of mass (CM) of the polymer coils derived from 
the CM-CM pair distribution function of the original sys- 
tem a new interest has emerged for the development 
of theoretical approaches able to predict these CM-CM 
correlations from first principles in a wide range of ther- 
modynamic conditions. Recently, a very simple route 
towards such theories has been opened with the deriva- 
tion, in the framework of the polymer reference inter- 
action site model (PRISM) Q, of an approximate rela- 
tion, involving only single coil structural quantities, be- 
tween the monomer-monomer and CM-CM pair correla- 
tion functions of a one-component homopolymer fluid Q . 
Indeed, with this formula, the predictions of which com- 
pare favorably with simulation data in wide density and 
solvent quality ranges [3|, |J, [fj( , any theory for monomer- 
monomer correlations can virtually be transformed into 
a theory for CM-CM correlations. This has been put 
to good use to derive effective CM-CM pair interactions 
starting from a very simple integral equation theory, with 
results in good agreement with simulation data, at least 
in the good solvent regime [ij . 

Lately, Guenza and coworkers have extended the re- 
lation derived in Ref. to binary polymer mixtures 
and, following a procedure analogous to that of Ref. [H, 
they have computed CM-CM pair distribution functions 
for melts and binary blends, with results in fair agree- 
ment with molecular dynamics data @ , thus demonstrat- 
ing the potential applicability of the theory proposed in 
Ref. [3[ to more complex macromolecular systems. In 
this Note, we seize the opportunity provided by this work 
and by the recent extension of the above coarse-graining 
strategy to diblock copolymer solutions^ , to present the 
generalization of the results of Refs. [3| and @ to arbi- 
trary mixtures of macromolecules of arbitrary structure. 

As in Refs. and the sought-for relation results 
from the application of the PRISM theory, now in its 
general multicomponcnt formulation, to macromolecular 
fluid systems in which two categories of sites are formally 
made to coexist: on the one hand, "real" or physical sites, 
which actually interact and whose interactions determine 



the microscopic structure of the fluid, and, on the other 
hand, auxiliary non-interacting sites which are only intro- 
duced to tag specific points of interest, typically following 
from a simple geometric definition like the CMs. 

In this approach, a macromolecular species x, present 
in the mixture at number density p x , is assumed to con- 
sist of n x classes of rigorously or approximately equiva- 
lent physical sites denoted by x = {ari, . . . , x„ x } and m x 
classes of rigorously or approximately equivalent auxil- 
iary sites denoted by X e {Xi = x nm+ i, . . . , X mx = 



}. We denote by N Xi , 1 < i < 



the 



number of sites of class Xi per macromolecule. 

Then, in reciprocal space at wave vector q, the PRISM 
equations read, in matrix form 



H(q)=W(q)C(q)[W(q)+H(q)]. 



(1) 



W(q) is the symmetric matrix of intramolecular struc- 
ture factors, or form factors, defined as 



wWq) = 



k=l 1 = 1 



(2) 



where Xi and Xj , 1 < i, j < n x + m x , are classes of real 
or auxiliary sites belonging to the molecules of species x 
and rj^ and r^ 3 are the positions of the sites of these 
classes on one macromolecule. W(q) is block diago- 
nal and, separating the coefficients according to the real 
or auxiliary nature of the involved sites, the diagonal 

blocks can be written ( ) . H(q) and C(q) 

are non-symmetric matrices simply related to the site- 
averaged total and direct correlation functions h XiVj (q) 
and c XiVj (q) through the relations 

II r ,, !q; = p x N x% N yj h x% y ] (a L ), C, iqi = p x c miVj (q), 

and, using the same separation of site classes as 
above, they are found to consist of blocks of the form 

/ H xy (q) H xY (q) \ d / C xy (q) C xY (q) \ reSDectively 
^Hxy(q) Hxv(q) ) VCx y (q) C XY (q)^' reS P eCT1Vel y- 

We now partially close the set of PRISM equations 
by setting all direct correlation functions involving auxil- 
iary sites identically to zero. This approximation is moti- 
vated on physical grounds by the fact that these sites are 
strictly non interacting and that all correlations in which 



2 



they are involved can only originate from intramolecu- 
lar correlations and direct interactions between physical 
sites. This results in a specific structure of the direct cor- 
relation matrix, consisting of simple blocks ( Cx y ( q ) ° ) . 

It is now a matter of simple matrix algebra to show 
that 

Hx y (q) = Wxx^W^H^q), (4) 
H XY (q) = Wxx(q)W^(q)H xY (q), (5) 

where W H '(q) denotes the inverse matrix of W xx (q), 
and, by transposition and multiplication by p y /p x , 



H yX (q) = H yx (q)W X3 J(q)W xX (q), 
H YX (q) = Hv^W-^Wxxfa). 

Combining these equations, it results that 



(6) 
(7) 



H XY (q) = W Xx (q)W- x 1 (q)H xy (q)W yy 1 (q)W yY (q). 

Note that there is an obvious simplification by p x on both 
sides of this equation, such that no explicit dependence 
in the densities is present. 

If the auxiliary sites are the CMs of monomer groups 
constitutive of the macromolecules, then Eq. JHJ is the 
generalized form of the expressions derived in Refs. @ 
and @. It is flexible enough to accommodate various 
ways of building coarse-grained representations of the 
macromolecules based on the CMs of their monomers. 
Indeed, for physical reasons or analytical convenience, 
one might choose to introduce the CMs of the sites of 
each class, leading to a coarse-grained representation 
with m x — n x sites for species x, as it was recently done 
for diblock copolymers (n x = 2) Q, or one might find 



useful to work with the global CM of the macromolecules 
only, hence m x = 1, or any intermediate grouping of non- 
equivalent site classes, then m x < n x . Of course, different 
choices for the various components of the mixture can be 
mixed. In the following, we explicitly discuss some of the 
simplest examples. 

The first one, relevant for homopolymer or colloid- 
homopolymer mixtures, corresponds to n x = m x = 1 
for all species. Then all matrices in Eq. ((8]) reduce to 
scalars and one finds 



W XlXl (q)W YlVl (q) 
W xlxl {q)W yiyi (q) 



N Xl N yi h XlVl (q), (9) 



as shown in Refs. i and 0, where a slightly different 
normalization was used for the form factors. 

More generally, for any species x comprising only one 
type of monomers (n x = 1) in an arbitrary mixture, 
any pair correlation function h XlV {(\) involving these 
monomers can be transformed into the corresponding 
pair correlation function hx iy {<\), involving the CMs of 
these monomers, using the simple relation 



ix iy 



(q) 



W XlXl (l) 



w x 



N Xl h Xiy (q), 



(10) 



' x lXl (q) 

where y denotes any type of site, real or auxiliary, on 
any macromolecular species, including sites of species x 
itself. 

We now consider the case of a one-component fluid 
with n x > 1 and m x — 1, i.e., the macromolecules 
are made of different types of monomers, but a coarse- 
grained description with their global CMs as the reference 
point is sought. Then, expanding Eq. ([8]), one finds 



J 



h Xl xM) = £ ^ix i (q)[W x ^]x i x 3 (q)^x i /ix,x fe (q)^jW- x 1 ] KfcK( (q)W' :c ,x 1 (q), 



(ii) 



r 



a result which might be useful to study copolymer or star 
polymer fluids. 

In conclusion, we have derived fully general approxi- 
mate relations between the monomer and CM pair corre- 
lation functions in macromolecular fluid systems. They 
extend previous results obtained for simple homopolymer 
fluids which have been found to compare very favorably 



with computer simulation data. With these equations, 
following the lines of the work done in Ref. [4| , a system- 
atic investigation of effective pair interactions in macro- 
molecular fluids seems achievable, at a modest computa- 
tional cost and with at least semi-quantitative accuracy. 
Work along this line is in progress. 
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